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Abstract 

The structure function of the vortex lattice of layered superconductor is calculated to one-loop 
order. Based on a phenomenological melting criterion concerning the Debye-Waller factor, we 
calculate the melting line of the vortex lattice, and compare our results to Monte Carlo simulation 
and experiment. We find that our results are quantitatively in good agreement with the Monte 
Carlo results. Moreover, our analytic calculation of the melting line of BSCCO fits the experiment 
reasonably well in a temperature range not far from T c . 



PACS numbers: 74.20.De, 74.40.+k, 74.25.Ha, 74.25.Dw 



I. INTRODUCTION 



Magnetic fields can penetrate into the layered high-T c superconductors(LHTS) and gen- 
erate the well-known vortex matter. 1 Thermal fluctuations and the effects of disorder are 
able to drive the vortex matter to undergo very complicated phase transitions among glass, 
solid and liquid phases.-^ This led to a burst of investigations both experimental and theo- 
retical, to understand the physical properties of such vortex matter.- One interesting aspect 
concerning vortex matter is to determine the phase transition line of the melting between 
the vortex solid state and the vortex liquid state.-^ Several significant experiments^^ have 
observed the phase transition of the melting between the vortex solid state and the vortex 
liquid state. In addition, magnetization jumps^*^ and specific heat spikes^ 1 ^ 1 ^ were ob- 
served, which indicate that the vortex lattice melting is a first-order phase transition. An 
often used theoretical description of the vortex lattice is the elastic theory,-^^^^ which is 
based on the lowest energy excitations on a perfect vortex lattice and can retain in fact most 
of the interesting physics. In the framework of elastic theory, one uses a phenomenological 
criterion, Lindemann criterion,— to calculate the phase transition line. 

Near T c vortices overlap and the elastic theory is questionable. Furthermore for vortex 
liquid, the elastic theory is not applicable. Another theoretical approach to study the phase 
transition in vortex matter is studying the thermal fluctuations of the more fundamental 
model, the Ginzburg-Landau(GL) model.— The model can describe very well both vortex 
lattice and vortex liquid in the region near T c . 

However, the full Ginzburg-Landau(GL) theory is too complicated, and one needs to use 
some sort of approximations to advance the theoretical investigation. Usually, the interested 
phase transition is located not far away from T c . Near T c , it is well known that one can 
use the lowest Landau level (LLL) approximation. 1,21 A number of researchers had studied 
the vortex liquid phase.— For solid phase, however, due to supersoft phonon modes, 
the perturbation theory was questioned. The problem was resolved by Rosenstein^ and it 
was found that all infrared divergencies are canceled. Non-perturbative Gaussian variational 
calculation had been carried out.— Spinodal line was determined and recently was confirmed 
by experiments.— 1 ^ By comparing the free energy of vortex liquid with vortex solid, the 
melting line had been obtained.-^ The result was used to explain the melting transition of 
YBCO type (not too high anisotropy) .—^ Recently it was also used to explain the melting 
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transition in low T c materials.— 

While for highly anisotropic ones such as the BSCCO, the case is far more complicated. 
The relevant model for highly anisotropic superconductors is the Lawrence-Doniach-GL 
model(LDGL).— . The LDGL model has been proposed to describe LHTS with weak Joseph- 
son interlayer coupling.-^^^ However the full studying of this model is not tractable. We 
assume that we can start from an effective LLL LDGL by integrating the higher Landau level 
modes and focus our studies on this model. The phase diagram of the LLL LDGL model 



had been investigated numerically in Ref. |3J, and the melting transition line was obtained. 



It is highly interesting if one can obtain the melting transition of the LLL LDGL model 
analytically. However even for the LLL LDGL, there is not a nonperturbative calculation 
of vortex liquid energy, therefore we have not yet determined the melting transition directly 
by comparing the liquid and solid free energy. 

Therefore we come back to use more phenomenological criterions to determine the melting 
transition. We might use the Lindemann criterion to study the melting transition, however 
unfortunately we do not know how to obtain the elastic modulii of the Lawrence- Doniach 
model and used them to obtain the melting transition. 

In this paper, instead we will use another criterion to determine the transition, the Debye- 
Waller factor criterion (DW criterion).—^ Based on this criterion, we will study analytically 
the model and the result will be compared quantitatively with the Monte Carlo(MC) one 
in Ref. 



34 and the model will also be used to calculate the melting transition line in BSCCO. 
The Debye- Waller factor is the number of the original height of the second Bragg peak 
(thermal fluctuation not considered) divided by the height of the second Bragg peak with 
thermal fluctuation considered. Due to thermal fluctuations, the Debye- Waller factor will 
be reduced. If the peak height is lowered to some extent, for example, 60%, the lattice 
will be melted. The applications of the DW criterion in both Yukawa system^ and three 
dimensional (3D) anisotropic case of the high-T c superconductors^ indicate that the criterion 
is quite accurate in determining the melting transition line. By using the DW criterion, 
we find that our calculations fit the result very well obtained by the MC simulation for 
studying the effective LLL LD model. The analytic calculation of the melting transition line 
is compared reasonable well with experiment.— We also stress that the calculations for the 
structure function and the Debye- Waller factor are fairly simple and straightforward. 

The paper is organized as follows. The model is described and a perturbative mean-field 
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solution is developed in section II. Then in section III the structure function of the vortex 
lattice is calculated to one-loop. A melting criterion is discussed in section IV, comparison 
with MC simulations and experiment was also discussed in this section. 



II. MODEL, MEAN FIELD SOLUTION, AND THE PERTURBATION THEORY 



A. Model 



We start from the following Lawrence-Doniach free energy: 
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where ipn is the order parameter defined on discretely labelled continuum layers, do is the 
layer thickness, d is the interlayer spacing, the covariant derivative is defined by D = V — 
i(27r/$ )> an d $o = (hc/2e). In the limit that do = d and d goes to 0, the LD model 
reduces to the 3D anisotropic Ginzburg-Landau model. For layered superconductors far 
from H c i (this is the range of interest in this paper), the magnetic field is homogeneous 
due to the overlap of the vortices. We choose the Landau gauge A = (By, 0,0), which 
describes a nonfluctuating constant magnetic field directed perpendicular to the layers. For 
simplicity, we assume a(T) = — a(l — t)T c ,t = T/T c , and other parameters are temperature 
independent. 

For convenience, within the LLL approximation, we use the following units to rescale the 
model:— the units of length of the "ab" plane and the "z" direction are = \Jh 2 / (2m ab aT c ) 



and £ c = a/ ft 2 / (2m c aT c ), respectively; the unit of magnetic field is H c2 , and the order 
parameter field is rescaled as ip 2 — > (2aT c /b')ip 2 . The dimensionless free energy in these 
units is 
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The dimensionless coefficient is uo = \f2Qm 2 t, where the Ginzburg number is defined by 
Gi = ^(32ire 2 K 2 ^T c 'y/c 2 h 2 ) 2 and 7 = a/ mj m a b is the anisotropy parameter. 



B. Mean field solution 



By minimizing Fql with respect to ip n , this standard variation problem leads to the 
well-known GL equation 



- a h 1p n + |^n|Vn = 0, 



(2) 



where = (1 — t — b)/2, b = B/H c2 , and TC = —(D 2 + b)/2. If is sufficiently small, 
the GL equation can be solved perturbatively. Within the LLL approximation, one gets the 
mean field solution of the GL equation 
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Where <£>(x) is the Abrikosov's lattice solution,— its definition is 
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and where a = \J 47r/v^3, = (|y?| 4 ) = J cell d 2 x\ip\ 4 (b/2n) rs 1.1596 is the Abrikosov's 
constant, the "cell" here is a primitive cell of the vortex lattice. Obviously, the mean field 
solution is independent of the layer index n. 



C. Fluctuation spectrum 

In order to get the excitation spectrum one expands the free energy functional around 
the mean field solution. The fluctuating order parameter ip n can be written as the sum of 
the mean field part and a small fluctuating part 

%l) n {x) = $(x) +Xn(x). (5) 

We emphasize here that the argument u x" in (5) stands for a 3D vector (i.e., x = (x, £3)), 
and the bold-face font (e.g., x) stands for the 2D vector in the "ab" plane in this paper. 
The field \ n can be expanded in a basis of quasimomentum eigenfunctions ip^ (within the 
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LLL approximation): 
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In order to do the perturbation calculation more conveniently (this can be seen later), we 
define (3-^ and 7k 



A< = (M 2 |y?k| 2 >, 

7k = ((¥»*) V-k^k), 



(7) 



while k = 0, (3q, 70 are shorted as f3, 7, respectively. We get 

y(x)-— [d*k e —^M (o, I 



(8) 



where fci, A: 2 G [—00, +00], &3 G [— ir/d, ir/d], d^ = exp [— and 7k = |7k| exp [i#k]- 
For simplicity, we have used in (8) the "real" field O k and A k , which satisfy the relations: 
01 = 0_fc, A\ = A_ k . Within the LLL approximation, at order ah, the eigenstates are Ok, 
A k , We find that it is convenient for us to get the eigenvalues by using <i k in the expansion 



of Xn(x). The eigenvalues are 



eo = eo + ^(1 - cos k 3 d) 

= a h ^-l + -/3 k + -|7k 
£a = £a + "4 (! - cos fad) 



- cos fc 3 d), 
' ^ In particular, 
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when A; — > 0, ~ 0.1a/i|/c| 4 , while eo has a finite gap. 



III. STRUCTURE FUNCTION OF THE VORTEX LATTICE 



In this section we calculate the structure function to order uo within the LLL approx- 
imation, i.e., neglecting higher correlations. Firstly, we calculate the density- density 
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correlation function defined by 



S(2,Z 3 ) = (p(x,X 3 )p(x+ Z,X 3 + Zg)) x , 



(10) 



where p(x) = \ip(x)\ 2 , and the subscript x here indicates the average over the unit cell. The 
correlation function is calculated using the well-known Wick expansion: 39 



Z 3 ) — S m { + CuSfluct, 



(11) 



where the first term is the mean field part, while the second term is the correction due to 
thermal fluctuations. 

A. Mean field contribution 



The mean field part is 



S m{ = (mx)\ 2 \<S>(x + z)\ 2 ), 



(12) 



The structure function is the fourier transform S*(q, 0) = J dzexp [iq ■ z]S(z, z 3 = 0), hence, 
the mean field part of the structure function is 



cfeexp [zq- z](|$(x)| 2 |$(x + z)\ 2 )^ 
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In order to derive Eq.(13), we have used the following relation: 



dx(/?(x)(p£(x) exp[— zx ■ q] 



4n 2 5 n (q — k) exp 



Til 



(n\ - rix) 



x exp 



q^ _ iq x q y ik x q y 
46 26 b 



(13) 



(14) 



where A is the sample area, and where we have used the notation: <5„(q) = J2 ni m 5(q — 
nxdx — n 2 d 2 ), rii = (l/27r)d! • q, n 2 — (l/27r)d 2 ■ q, di, d 2 are the reciprocal lattice basis 
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vectors 

which are dual to the lattice basis vectors 

di = (a/ Vb, o) ; d 2 = (a/2 A . (16) 

B. Fluctuation contribution 

We calculate the fluctuation contribution of the structure function (S COTT .) to one loop. 
For convenience, the results are divided into four parts: 

5 cor ,(q,0) = 5 1 (q,0) + 5 2 (q,0) 

+53(q,0) + 54(q,0), (17) 

where Si(q, 0) is the fourier transform of {$(x)$(x + z)Xn( x )Xn( x + z ) + c.c.) x , S^q, 0) 
is the fourier transform of ($(x)$*(x + z)Xn( x )Xn( x + z) + c.c.) x , S^q, 0) is the fourier 
transform of (\&(x)\ 2 \x n ( x + z)\ 2 + |<&(x + 2;)| 2 |x„(a;)| 2 ) x , and S^q, 0) is the fourier trans- 
form of j^-(\ip(x)\ 2 \(p(y)\ 2 ) x (is 2 ). We emphasize that the final term is due to the vacuum 
renormalization, which cause the shift v in ip n (x) = vip n (x) + Xn( x ) be renormalized. To 
one-loop order, let v 2 = z/ 2 + uu 2 , here v\ = a^j Pa, and v\ is given by minimizing the 
effective one-loop free energy 
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where L x , L y , L z are the scales of the sample. From Eq.(18), we get 
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Each term of the r.h.s of (17) is given as follows: 
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where q = k + Q, k is the fractional part of q, while Q is the integer part of q. After 
combining the mean field part and the one loop correction part of the structure function, 
we get 
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It is very interesting to notice that each of the four terms Si(i = l,---,4) is divergent, 
respectively, as k — > 0, however, the sums Si, S 2 and 53, S 4 are not. Here we just take the 
sum Si + 62 as an example: 



Si(q,0) + S 2 (q,0) 
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exp 
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In order to see it more clearly, we use \/b to rescale the momentum. As it can be shown 
that k x k y + #k = 0(k i ) when k — > 0, the function (A^/fy + (k x Q) z + # k ) — > (k x Q) z , 
and 1 — cos(k x k y + (k x Q) z + # k ) -»(kx Q)^, hence it will cancel the 1/k 2 singularity of 



^2/{l A + d 2 l 2 A /2). It is also easy to show that S3 + S 4 is the case. Consequently, we get 
a not divergent result. The fluctuation correction of the structure function (for non-peak 
region) is shown in Fig.l. 
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IV. MELTING OF THE VORTEX LATTICE 



A. A melting criterion 

The above calculations also indicate that thermal fluctuations will reduce the intensity 
of the Bragg peak of the structure function. In fact, the Debye- Waller factor has been 
used to describe the melting of the lattice system.-- If the intensity of the peak is lowered 
to some extent, for example, 60%, the lattice will be melted. In 3D case, we know the 
exact melting transition temperature via a different method. 2 - With this temperature, we 
find that the one loop calculation of the Debye- Waller factor is reduced to 50%. This does 
not mean that the criterion "60%" is wrong as the higher order correction to the one loop 
calculation usually will increase this value to some number above 50% (it is too complicated 
to calculate the structure function to two loop and we will leave it to future studies). In 
this paper, we also calculate the structure function to one loop order. Thus we will use the 
"one loop criterion" that the Debye- Waller factor calculated to one loop is reduced to 50% 
at the melting transition line. 

According to the definition of the Debye- Waller factor, we denote the ratio of the one-loop 
value to the mean field value of the intensity of the second Bragg peak by p, we get 

9 ~ ' (29) 

where Qi denotes the shortest reciprocal vector of the triangular vortex lattice. We define 
the critical value of p corresponding to melting by p c . According to the above discussion, 
the one loop criterion corresponding to p c is about 50%. 



B. Comparison with MC simulations 



Now we c omp are our results with MC simulations of the LLL layered system in Ref. 



34. In Ref. 



34 



, the authors use two dimensionless parameters g and rj to describe the 



melting transition of the system. The g and 77 measure the intralayer and interlayer coupling, 
respectively, and they are equivalent to the t and b in this paper (in fact, they are the 
functions of t and b). In order to carry out the comparison, first, we make our notations 



consistent with the ones of Ref. 
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9 = aB ^m^B = ~ ah y^ (31) 

h 2 l 

V = 2m c {di c ) 2 \a B \ = 2d 2 a h (32) 

Furthermore, from Eq.(26), (27), it is easy to see that /2(Qi) + fz is only dependent on 
d 2 cih, we define f(d 2 dh) = /2(Qi) + /3, then, we have 

f(dW) = ( P -l)^%. (33) 

After combining Eq.(31),(32),(33) together, we get 

1 „, 1 . , , 87iq 2 , , 

-w/^-^-^it- (34) 



According to the knowledge from 3D anisotropic model,— we choose p c = 0.475, hence, 
Eq.(34) gives out the relationship between g and rj. 



In Ref. 



34j, the MC simulation was employed to determine the melting transition. On 



34 



the melting transition the value of g is a function on rj. We denote the function g in Ref. 
by gHM, and the g determined by Eq.(34) is denoted by g. The result of the comparison is 
given in Table I. For typical LHTS such as BSCCO, £ a b is about 25A, 7 is about 200, H c2 
is about 50T, T c is about 90K, and the interlayer spacing d is about 4A. For temperature 
and magnetic field at 75K and 400G respectively, the value of 77 is about 0.01 (actually 77 
on the all points on the theoretical curve on fig. 2 is less than 0.1). when 77 increases, \g\ 
is gradually larger than \guM\- As discussed in Ref. |34j that the finite size effects of MC 
simulation become stronger as rj increased, and the finite size effects lead the |<7hm| to be 
less than its actual value. In summary, we find the two results fit very well for not too big 
77 (less than 0.1) . This also demonstrates that the DW criterion works well. 



C. Comparison with experiment 



In Ref. 



371 . the material parameters describe BSCCO: k = 100, 7 = 270, T c = 86K, 
d = 15A. The interlayer spacing do and H c2 have not given, we find d = 4.1 A, H c2 = 50T 
give the best fit to the experimental data, the results is shown in Fig. 2. The deviation become 
large as T reduces, this is expected as the effects of disorder will be enhanced as temperature 
lowed. The effects of disorder tend to lower the curve. However, the effects of thermal 
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TABLE I: Comparison of guM with g 



7] 


guu 


9 


n 005 


-5 9 


-4 75 


0.01 


-4.5 


-4.36 


0.02 


-3.9 


-3.96 


0.06 


-2.98 


-3.31 


0.10 


-2.74 


-3.02 



fluctuations dominated in the region of our interest (near T c ). The comparison indicate that 
the effective LLL LD model is quite good to describe the melting phase transition of the 
LHTS near T c . In future work we will include the disorder effect, and we expect that the 
result can be extended to the region with lower temperature. 

V. CONCLUSIONS 

To conclude, we have calculated the structure function of layered superconductors and 
the melting line can be obtained quantitatively by "one-loop" DW criterion,i.e. the ratio of 
the one-loop value of the intensity of the second Bragg peak of the structure function to the 
mean-field value is about 50%, the solid melts. With this criterion, we calculate the melting 
line and compare the results with existing MC results and experiments. Our results fit the 
MC results very well. Moreover, our results fit the experimental data reasonably well in the 
range not far from T c (for BSCCO, T c = 86K, the range we fitted is from 72K to 86K). 
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Figure Captions 

Figure 1 Fluctuation correlation to the structure function of the Abrikosov vortex lattice. 
The peaks at reciprocal lattice points are removed, only the correction to the non-peak region 
is plotted (i.e. only /i(q) is plotted. 

Figure 2 

Comparison of the theoretical melting curve(line) of highly over doped BSCCO with 
experimental results. 
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